Ab-initio calculations of the optical properties of the Si(113) 3x2ADI surface 
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We investigated the stable silicon (113) surface with a 3x2ADI reconstruction by ab-initio meth- 
ods. The ground state properties have been obtained using the density-functional theory. We 
present the dispersion of the electronic band structure, where the surface bands have been distin- 
guished from the projected bulk bands by calculating their localization in the slab. The optical 
spectra, here the reflectance anisotropy (RAS), have been obtained within the independent particle 
random phase approximation. We identifled surface features in the spectra tracing them back to the 
responsible electronic states and, studied their localization in the slab. A comparison with available 
experimental data for the band structure and the RAS shows a good agreement. 

PACS numbers: 78.68.-l-m 68.35.Bs 73.20.At 
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I. INTRODUCTION 

Although being of high index, the vicinal Si(113) sur- 
face is one of the stable silicon (Si) surfaces.— This sur- 
face is of technological interest, since atomically smooth, 
ultrathin oxide films can be grown on iti^ Hence, it is 
dealt as a candidate for the next generation of waferSjS, 
and has potential applications in nanostructures tech- 
nology. Si(113) can be used as a substrate for the self- 
assembled growth of germanium (Ge) nanodots^ and 
nanowire9^i2', as well as Ge^iS and SiGe^° islands. Be- 
sides the technological interest, Si(113) is also of funda- 
mental one, because it shows phase transitions between 
its 3x2 and 3x1 reconstructions. These phase transitions 
can be induced by temperatur o^^i^^'^'^:^^'^^:^^ and by 
contamination,iii^iiSi2ai21i22i2^i2i The latter ones affect 
also the discussion about the surface reconstruction. In 
spite of the measurement of a 3x 1 reconstruction,^^ most 
of the experiments show at room temperature (RT) a 
3x2 surface periodicity^^'ii^'^o^'Sl^^ The finding 
of the 3x1 periodicity might be due to contaminationii 
and the 3x2 one is assumed as the surface unit of clean 
Si(113) at RT. 

Various surface reconstructions have been proposed 
for the Si(113)3x2 surface. The most probable sur- 
face model is the ADI (adatom-dimer-interstitial) re- 
construction of Dabrowski,— which is in agreement 
with experimental high-resolution transmission electron 
microscopy (HRTEM) ^^'^^ and scanning tunneling mi- 



croscopy (STM) 
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images. It has been con- 



firmed by other theoretical investigations, tooi^i^ii^i^^ 
A top view of the surface unit cell is shown in Fig. [T] 
The surface reconstruction shows two adatoms, and two 
pentagons with a dimer along the [110] direction, with 
one pentagon hosting an interstitial atom in its center. 

Earlier structured models were based on the structure 
proposed by Ranke^^ Si(113)3xlAD, which is similar to 
the ADI one. It can be deduced from the ADI model by 
eliminating the interstitial atom and relaxing the atoms 
in z direction. Despite being intrinsically 3x1, also in 
the 3xlAD model the 3x2 periodicity can be obtained 



either by canceling one pentagon in each 3x2 surface unit 
cell ( "surface- void" ^33] ) , or by manipulating the z coor- 
dinates of the pe ntag ons in order to break the symme- 
try ("corrugated" [If, "puckered" [Hil , "shifted" [H, 
etc.). Since the 3x1 periodicity is also interesting for 
the phase transition, being the periodicity of the high- 
temperature phase, besides the 3xlAD model also a 
3xlAI one, with an interstitial atom at each pentagon 
has been proposedi^^i^ However, also other surface mod- 
els have been assumed, generally based on the models of 
Ranke or D§,browski42ji^ 

From the theoretical point of view, the most stable sur- 
face (i.e., the one with the lowest surface energy) should 
be selected. For some of the models above, the sur- 
face energy has been calculated, unfortunately by using 
different methods. There exist some studies assuming 
more than one surface reconstruction^SiiHi^iiid^ within 
the same framework, but none of them considered all 
available surface models for a comparison. Most of the 




FIG. 1: (color online) Top view_ of the Si(113)3x2ADI ceU 
(a::, y plane corresponding to [lIO] and [332], respectively). 
Silicon atoms are in large dark spheres (red). The interstitial 
Si atom (atom 18) is depicted as a large light (yellow) sphere. 
Hydrogen atoms are small light spheres (blue). For a side 
view see Figs. [S] and [S] Surface atoms (see sketch) include 
two pentagons with dimers (atom pairs 1-2 and 5-6) and two 
adatoms (atoms 3 and 9). 
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calculations have shown that the Si(113)3x2ADI surface 
has the lowest surface energy with respect to the other 
assumed reconstructions. 

Nevertheless, not all measurements agree with the ADI 
model.— i^!^ The main criticism to the ADI model is 
due to the presence of a mirror plane, which should 
yield the extinction of some diffraction reflections,-^ be- 
cause the corresponding patterns have been seen in some 
experiment.— Thus, models different from the ADI one 
are also being considered and a final conclusion has not 
been drawn yet. 

An alternative approach for surface investigations is 
given by the study of its optical properties. A very effi- 
cient technique is the reflectance anisotropy spectroscopy 
(RAS), which can be used in situ for the surface determi- 
nation. The interpretation of such spectra is a demanding 
task. For this purpose it is necessary to have theoretical 
reference data together with a correlation of the spectral 
features and the surface reconstruction elements. 

The RAS of the Si(113) surface has been measured 



since ten years 
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The first theoretical investigations 



of the RAS for the Si(113) surface has been performed in 
1998 by using a dipole approximation for the ideal bulk 
atomic positions.-^ In the same year, a tight-binding cal- 
culation has been done for the Si(113)3x2ADI model^^, 
where the bulk derivative- like peaks have been found with 
the wrong sign. Although tight-binding calculations yield 
good results for bulk silicon, the low-energy range of the 
spectra, which is related to the surface, is usually not 
described very well.^^ 

For this purpose, we choose an ab-initio approach to 
investigate the optical properties of the Si(113) surface. 
The goal of this study is to find out whether surface- 
related structures in the spectra can be traced back to 
some characteristic surface reconstruction elements of 
Si(113)3x2ADI (pentagons, interstitial, adatoms, etc.), 
and which kind of changes can be predicted with respect 
to other surface models. In this way optical spectroscopy 
could be used to discriminate between different models 
for reconstruction of the Si(113) surface. Due to the wide 
acceptance of the Si(f 13)3x2ADI, we used this model as 
a starting point. 

This article is organized as follows: After summarizing 
shortly the methods used for the calculation (Section lTT)) . 
we first focus on the ground state properties (SectionlIII|). 
Since the optical spectra are determined by the electronic 
transitions, an investigation of the electronic structure is 
presented in Section IIVI A detailed discussion of the 
optical spectra is given in Section |Vl and, finally, we 
summarize and draw a conclusion. 



II. METHOD 

We have performed ab-initio total energy calcula- 
tions using the periodic slab method within the frame- 
work of the density-functional theory^^ (DFT) as im- 
plemented in the ABINIT^^ a,nd TOSCA^i packages. 



For the exchange-correlation energy in the Kohn-Sham 
equations^^ the local-density approximationSSi-^^ (LDA) 
has been chosen. The eigenfunctions have been expanded 
into plane waves using pseudopotentials, here normcon- 
serving ones in the TrouUier-Martins style^^ 

After converging the ground state structure with 
ABINIT, the optical properties have been calculated us- 
ing the TOSCA packaged The probability P^.ck of the 
electronic transitions between the valence (v) and the 
conduction (c) states with electronic eigenenergies Ecu 
and iJ^k for light polarized in j (j = x, y, z) direction at 
a given point k in the reciprocal space has been calcu- 
lated as the diagonal elements of the velocity operator.— 
Working within the independent particle random phase 
approximation^ (IPRPA), local-field, self-energy, and 
excitonic effects are neglected. For our system we assume 
that they can be well described by the scissor operator 
approach;^ hence we stick to the IPRPA to describe the 
optical properties. In order to obtain the RAS, in a first 
step the imaginary part of the slab polarizability has 
been calculated by 
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for an energy Here, A is the surface area, m and 
e the electron mass and charge, respectively. In order 
to smooth the resulting spectra, a Gaussian broadening 
has been applied, because it leads to more well-defined 
structured spectra than the Lorentzian one. The RAS is 
defined as the difference of the deviation from the Fres- 
nel refiectivity ARj/R for the orthogonal polarizations 
within the surface plane (here the x and y direction): 
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The deviation from the Fresnel refiectivity for normal 
incident light is given by^ 



AR^ 
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where ajj'^^ is the complex surface polarizability and 
ab(w) the bulk one (cq: velocity of hght). The surface 
polarizability can be obtained by subtracting the bulk 
contribution from the slab one. Since the imaginary part 
is taken from the ratio, the bulk part in the numerator 
vanishes, and Eq. ^ holds also for a^^^ instead of ct^j'^^ ■ 
The complex polarizability function a is derived from the 
imaginary part via the Kramers-Kronig transform. 

Finally, to keep into account the self-energy and exci- 
tonic effects in an approximative way, a scissor operator 
shift^^ has been applied to the eigenenergies. 
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III. GROUND STATE 

For our calculations we have chosen the cell with the 
3x2ADI surface reconstruction as described above. We 
have used slabs with 11 and 7 double layers (DL) of Si, 
respectively. The amount of vacuum was taken equal to 
the slab thickness. Each DL consists of 12 Si atoms, just 
the surface DL contains one atom less. Since the slab is 
not symmetric, the dangling bonds at the bottom of the 
slab have been saturated with hydrogen (H). Therefore, 
the cells contain 131 or 83 Si atoms, respectively, and 18 
H atoms. A top view of the slab is shown in Figlll while 
the side views of the 11 DL cell can be seen in Figs. [5] 
andm 

The convergence of the total energy requires 4 k 
Monkhorst-Pack^^ points in the irreducible wedge of the 
Brillouin zone (IBZ), and a kinetic-energy cutoff of 16 Ry. 
With these parameters, the error in the total energy is 
less than 0.046 eV per atom. 

A relaxation of the topmost 6 DL, computing the forces 
acting on the atoms (the remaining DL have been kept 
fixed to bulk positions) shows significant changes with 
respect to the bulk positions just in the topmost 4 DL. 
Therefore, these 4 DL have been optimized till the forces 
are less than 0.08 eV/A. 

We have calculated the surface energy by subtract- 
ing the energy of the the hydrogenated surface and the 
bulk energy from the energy of the slab. The result- 
ing value of 12.74 eV per unit cell agrees very well with 
the value of 12.69 eV obtained by Stekolnikov et el4^ 
As a consequence, the value of surface energy per unit 
area of 87.87 meV/A^ is in good agreement with the one 
of Stekolnikov et al.-^- (87.36 meV/A^). However, it is 
slightly lower than the ones obtained from other groups 
(97 meY/A^fsd^, 90.4 meV/A^fsl], 94.8 meV/A^gl) 
using various methods. Nevertheless, it is lower than that 
for other surface reconstructions and the overall agree- 
ment is good. 

IV. ELECTRONIC STRUCTURE 

After the ground state has been converged, a first step 
towards the optical properties is an analysis of the elec- 
tronic structure. We started with the calculation of the 
electronic density of states (DOS) as displayed in Fig. [21 
where the highest occupied band (HOMO) is set to eV. 
There, a comparison of the DOS calculated with the 7 DL 
and the 11 DL slab is shown. The DOS has been obtained 
using a Brillouin-zone summation over 25 k points, to- 
gether with a Gaussian broadening of 0.06 eV. The gen- 
eral shape is similar to the one of bulk silicon for the 
occupied bands. In addition, around the energy gap 
surface-related structures are visible: one strong peak 
at the conduction band edge and a peak with a shoul- 
der and another shoulder close to the bulk-like part of 
the DOS at the valence band edge. The use of the 7 DL 
and the 11 DL supercell leads to just minor differences. 
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FIG. 2: (Color online) Electronic density of states for the 
11 DL cell (solid dark-green lines) and the 7 DL cell (dashed 
light-blue lines). The zero of the energy scale is set to the 
top-valence state, and the height of the DOS is scaled with 
respect to the total number of atoms in the corresponding cell 
(see text). 

By normalizing the DOS to the total number of atoms 
in the slab, the height of the surface peaks (around the 
gap) appears to be smaller for the thicker slab. 

The valence peaks around the Fermi energy, which 
are in addition to the bulk DOS, have been also found 
experimentally and they correspond to surface bands. 
Using angle-resolved ultraviolet photoemission spectra 
(ARUPS) Myler and Jacobi™ found two peaks sepa- 
rated by 1.6 eV, where the second one is broad and 
could hide two single peaks. A similar experiment has 
been performed using angle-resolved photoelectron spec- 
troscopy (ARPES) showing two single peaks separated 
by 0.4 eV, and a third one at lower energy with a non- 
dispersive character, which is not resolved in the other 
experiment F^^ii^i^ At normal emission, the two close 
peaks fall together resulting in a single broad peak, which 
has been found also by Myler and Jacobi."'^^ Thus, the sur- 
face peak around the HOMO in the DOS is well defined 
in the experiments, and it contains probably two peaks 
as indicated by the shoulder. Furthermore, the core level 
spectra of Hwang et alJ^^ were fitted using three sur- 
face peaks, where the second shoulder in our DOS might 
refer to that third peak. In scanning tunneling spectra 
(STS)^ two peaks are visible. The low energy maximum 
has been attributed to the tetramers and the high energy 
maximum to the adatoms of the Si(113)3x2ADI recon- 
structed surface. For the peak at the lowest unoccupied 
band (LUMO), which should be accessible by inverse pho- 
toelectron spectroscopy, no experimental results exist. 

After identifying the surface structures in the DOS we 
have investigated the dispersion of the electronic eigenen- 
ergies using the 11 DL supercell. The results are dis- 
played in Fig. O In order to separate the surface bands 
from the bulk background we have calculated the local- 
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FIG. 3: (Color online) Dispersion of the electronic eigenen- 
ergies along the high-symmetry directions of the irreducible 
part of the Brillouin zone (IBZ). According to the localization 
of the corresponding states in the slab, light dots correspond 
to bulk states and dark ones to surface states. The dot size 
accords to the value obtained by Eq. Q (see text). 

ization L„ for each state n by 

L„ = y drrnkir)0{z)i^ni. (4) 

where 0{z) is a realspace boxcar function^^, meaning 
6{z) — 1 for a given range in z direction and zero oth- 
erwise, ipnk is the single band wave function for a given 
point k and the integral is taken over the whole super- 
cell. Considering the range, e.g., from the subsurface DL 
to the vacuum above, the value of L„ corresponds to the 
integrated charge for that state located in the x-y slice 
with the chosen z range. In this way, this value indicates 
the amount of localization in the range selected by 9{z). 

Fig. [3] shows the results. The surface-localized bands 
are drawn with dark symbols (the z range has been cho- 
sen to cover the topmost two DLs), and the bulk- localized 
bands with light symbols (the z range has been chosen to 
cover the lowest seven DLs where the relaxation has not 
changed the bulk positions). The size of the dots corre- 
sponds to the amount of localization in this area. This 
means, e.g., for the topmost four valence bands at the T 
point, that the top band is less localized in the bulk area 
than the second one and the third band is less localized 
in the surface area than the fourth one. 

With this procedure we have identified three surface 
valence bands (So, Si, S2) and two surface conduction 
bands (Ci, C2), as visible in Fig. [31 All surface bands 
show a partial overlap with the bulk band structure. 
This points to a reduced influence of excitonic effects, 
and thus, the IPRPA scissor approximation is expected 
to work well. The state S2 appears as a state folded at 
J along the [110] direction, and hence shows 3x1 peri- 
odicity. The other surface states show the 3x2 one. Of 
course, the eigenvectors have 3x2 symmetry for all sur- 
face states due to the choice of the 3x2 supercell. In both 



the experiments, ARUPS and ARPES, valence surface 
states with 3x1 and 3x2 periodicity have been found, 
but with a larger energy difference, which is here between 
J and K at about 0.16-0.19 eV. The shape of the bands 
Si and S2 are in excellent agreement with the measured 
ones, and show a similar dispersion of about 0.17 eV and 
0.20 eV, respectively, compared with the experiments of 
An et al,^ (Si: 0.15 eV and S2: 0.30 eV) and Myler and 
Jacobii^ (Si: 0.20 eV and S2: 0.15 eV). Hence we con- 
clude that these two states are the ones appearing in the 
experiments. Along the F-J direction the separation of 
the bands shows large variations within the two experi- 
ments, where we have here an overlap. The discrepancy 
might be due to the negligence of quasiparticle correc- 
tions as well as due to experimental setups. Consider- 
ing available theoretical investigations of the dispersion 
of the eigenenergies of the Si(113) surface (besides the 
tight-binding study of Wang et alj^ using a puckered 
surface model, which is hence not comparable), there ex- 
ist just the investigation of Stekolnikov et ali^. Within 
their results they obtain also the surface states So and 
Si, where the S2 might be hidden in the bulk-projected 
band structure. The energy spacing between these states 
is similar to the one obtained here. 

By calculating the squared modulus of the wave func- 
tions of the surface states we have been able to local- 
ize them more precisely in the slab (see Fig. 2]). The 
states Sq and Si are located at the adatoms (atom 3 and 
atom 9, respectively, see Fig. [T|) and show a dangling 
bond character, which is in agreement with the results 
of Ref. [3^ and the assumption of Ref. . The states 
Ci and C2 have been found at the atoms 4-7-8 and 5-6- 
7-8 of the pentagon with the interstitial, which is stated 
similarly in Ref. [s^. Furthermore, the surface state S2 
is located at one adatom and two atoms of the empty 
pentagon (atoms 1-2 and 3), which has not been inves- 




FIG. 4: (Color online): Top view {x-y side like Fig. [T| of 
the Si(113)3x2ADI cell together with the charge density of 
the surface states So, Si, S2, Ci, and C2 (see text). Note 
the overlap between the localization of some of the surfaces 
states. 
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FIG. 5: (Color online): Side view {x-z plane) of the 
Si(113)3x2ADI cell together with the charge density of the 
low-energy Tamm-like state. Si atoms are drawn as large 
spheres, H atoms are as small spheres. 



tigated before. However, this was already assumed by 

Arabczyk et al.— . Anyway, also this state is present in 
experimental^'^'21i30,33^ theoretical22i33^i36 gr^^ 

images. 

Besides the surface states above, we have obtained an 
additional valence surface state at very low energy cor- 
responding to the separated peak at around —12.5 eV 
in the DOS of Fig. [21 This state shows a dispersion of 
less than 16 meV, which follows the dispersion of the 
lowest valence states, which is 20 times larger. Thus, 
this nearly flat state can be considered as a Tamm-like 
state.-Si This state is located in between the atoms of the 
pentagon with interstitial (atoms 4-5-6-7-8 and 18), as 
displayed in Fig. [5l The large separation between these 
sites in the real space is responsible for the flat character 
of the band. However, it is not the first time that such a 
state has been found. Pandey et al.^^ found a flat surface 
band at low energy (« 10.8 eV) which is partially below 
the bulk valence bands for Si(lll):SiH3. Also a theoreti- 
cal and experimental investigation of the GaAs(110)lxl 
surface^ shows, that one surface state is located below 
the valence band minimum with a gap in between. The 
Tamm-like state found in our calculation should be ac- 
cessible within measurements and its experimental find- 
ing would confirm the validity of the surface model used 
here. Due to its low-energy character, transitions from 
that state will not appear in the low or middle energy 
range of the RAS. 



V. OPTICAL PROPERTIES 

In this section we present results for the optical spec- 
tra, in particular for the RAS, of the Si(113) surface. We 
focus on the convergence with respect to the k points 



and the slab thickness, we analyze the RAS regarding 
surface contribution and responsible transitions, and we 
compare with available experimental results. 



A. Convergence tests 

For the optical properties we have performed conver- 
gence tests concerning the number of conduction states 
and the number of k points used in Eq. ([1]). A number 
of 460 bands in total (189 conduction states) have been 
sufficient for the 11 DL cell as well as 360 bands (185 
conduction states) for the 7 DL one, respectively. Each 
state has been expanded into plane waves. The use of 
4000 G vectors over the available set of Ri24000 plane 
waves has been enough for the convergence of the ma- 
trix elements. The convergence tests with respect to the 
special k points are more demanding. In our case, the k 
points can be chosen in one quarter of the BZ correspond- 
ing to the IBZ. A Gaussian broadening of 0.08 eV gives 
reasonably structured spectra. We have performed tests 
with various sets of k points for the 7 DL cell as shown 
in Fig. [S] As visible in the figure, up to 25 k points 
the spectra still show some deviations. With a single k 
analysis^^ it has been figured out that the convergence 
in the x direction is slower than in y direction and thus, 
an increase of the number of k points just in x direc- 
tion (resulting in the set of 55 k points) yields a strong 
improvement of the convergence. Comparing the spec- 
tra based on 25 and 55 k points we have concluded that 
convergence has been already achieved with 25 k points. 
Incidentally, using the F point only, a strong anisotropy 
appears (see inset of Fig. ^ which can be traced back to 
infinite chains in x direction of the bulk layersi^. 

Since the slab thickness can influence the spectra, we 
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FIG. 6: (Color online): Convergence of the RAS of 
Si(113)3x2ADI with respect to the k point summations. The 
1 k point sets (F and Baldereschi) are shown in the inset, since 
the corresponding RAS is larger by one order of magnitude. 
Convergence requires 25 k points in the IBZ. 
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FIG. 7: (Color online): Convergence of the RAS of 
Si(113)3x2ADI with respect to the slab thickness (7 DL, 
11 DL, 15 DL) using only the Baldereschi point in Eq. ([T]). 
Convergence is satisfactory using the 11 DL slab. 



have calculated the RAS using only the Baldereschi^^ 
point for various cell sizes, i.e. using slabs of 7 DL, 11 DL, 
and 15 DL. The 15 DL cell has been created by adding 
four bulk DL to the 11 DL slab. The outcome of the 
test is displayed in Fig. [71 It is clearly recognizable that 
the 7 DL slab does not describe the low energy region 
of the RAS correctly, where the differences between the 
11 DL and the 15 DL cell are just minor ones. Hence, 
the use of the 11 DL cell produce reliable results for the 
surface spectra. Note that the spectra using only the 
Baldereschi point are not converged, but the use of just 
one point is sufficient for the check performed here. It 
has been already shown that the convergence regarding 
the thickness improves using more k points.— Concerning 
the k-point convergence, we expect that our 25 k-points 
set, which has been shown to be enough for the 7 DL 
slab, will also be at convergence for the 11 DL one. 



B. Analysis of the RAS 

We proceed now to analyze our fully converged RAS 
spectra (using the IIDL supercell and 25 k points in the 
summation) with respect to the surface structure. Thus, 
we have performed a layer-by-layer decomposition as de- 
scribed in Refs. [bgIIt^It^ I. For this purpose, a real space 
cutoff has been used in order to separate the contribu- 
tions coming from defined layers of the slab by intro- 
ducing a boxcar function 0(z) in the calculation of the 
dipole matrixelements. The modified transition proba- 
bility PJij. ^1^. can be rewritten as 

PL,c^--^^ J dr^:k^(^)^V'ck (5) 



and the summation of Eq. ([T]) changes to 

xd{Eck ~ Eyk - huj) 

Note that only one of P and P* must contain the 6 func- 
tion. The decomposed spectra are presented in Fig. [10] 
together with a side view of the slab. We have num- 
bered the DL from the top to the bottom. For the anal- 
ysis we have chosen three regions: The topmost two DL 
(DLOl and DL02 (surface), while DLOO is the vacuum 
layer above), the third and the fourth DL (subsurface, 
DL03 and DL04) where the optimization changed the 
atomic positions with respect to the bulk ones, and the 
bulk part together with the hydrogens (bulk, DL05 till 
DLll, while DL12 is the vacuum layer below), since the 
relaxation did not affect the atomic positions in DL05 
and DL06. 

In the low energy range of the spectra various peaks 
are visible. For further discussion we have labeled the low 
energy peaks with PI, P2, P3, P4, and P5 (see Fig. [TU)) . 
Looking at the decomposed spectra in Fig. [8| one can see 
that the the peaks in the full spectra (containing the con- 
tribution of all DL) up to 2.5 eV (P1-P5) can be traced 
back to the surface and the subsurface DL, where the 
structures at higher energies are also bulk determined. 
Nevertheless, also at energies larger than 3 eV there ex- 
ist surface related structures. However, these structures 
are due to surface-bulk, bulk-surface transitions and sur- 
face resonances. PI and P3 originate from the surface 
layers only, whereas the peaks P2, P4, and P5 contain 
contributions from both the surface and the subsurface 
layers. Thus we can conclude that for an analysis of PI 
and P3 only transitions between surface bands need to be 
considered, whereas for P2, P4, and P5 also resonances 
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FIG. 8: (Color online): Side view of the 11 DL slab (left side) 
and layer- by- layer analysis of the RAS (right side). Besides 
the decomposed spectra (solid lines) also the full spectra (dot- 
ted lines in the upper three panels and solid line in the lower 
one on the right ) is drawn for comparison. The arrows indi- 
cate where the contributions to the spectra come from. The 
DL assignment is explained in the text. 
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FIG. 9: (Color online): bulk coefficients A and B from Eq. ([7|. 
For energies smaller than 2.5 eV there is no contribution of 
Re(47ra') to the RAS. 



or high-energy bulk states have to be investigated. 

After identifying the surface peaks in the spectra, we 
want to identify the responsible transitions. Since for 
the deviation from the Fresnel reflectivity Eq. Q the 
whole complex polarizability is considered, the RAS is 
also influenced by the real part of a. Thus, a overlap of 
spectral structures coming from the imaginary and the 
real part is possible. A feasible procedure for determining 
this overlap has been described in Ref . [ill , where Eq. ^ 
was rewritten as 

A R 

= A Im[47ra^-j] + B Re[47ra5^] . (7) 

The energy-dependent coefficients A and B are deter- 
mined only by the dielectric function of the bulk crystal. 
With the magnitude of A and B one can figure out if the 
RAS is mainly determined by Im(47rQ;*') or by Re(47rQ;*'). 
In this way, A and B are surface independent coefficients. 

Because of this, we analyzed these coefficients for our 
system, in detail, we determined A and B for bulk silicon. 
They are universal for all silicon surfaces, independent of 
the orientation and reconstruction. The result is shown 
in Fig. [HI Besides A and B, also \A\ is plotted in or- 
der to compare the magnitude of the coefficients. The 
result is the following: Every RAS of any clean silicon 
surface is determined only by Im(47ra^) for an energy 
up to 2.5 eV and by both, Im(47ra^) and Re(47rQ:'^), for 
higher energies. In detail, it is mainly determined by 
Im(47ra") for 2.5-3.2 eV, by both equally from 3.2-3.4 eV, 
mainly by Re(47rQ;'') for 3.4-5.1 eV, and by — Im(47ra'^) 
for higher energies. Thus, the surface peaks P1-P5 can 
be traced back to the differences just between Im(47rQ;|^) 
and lm{4:Trayy) . These are the differences we investigate 
now in order to find the origin of the surface features in 
the RAS. 

Since the RAS in the low-energy range (up to 1.3 eV, 
peaks PI to P3) has been reproduced already sufficiently 



well using the Baldereschi point only, based on this spec- 
tra the contribution of various transitions to the first 
three surface peaks has been analyzed. For the other 
spectral structures all k points have been taken into ac- 
count. 

The first peak PI of the RAS displayed in Fig. [TO] ap- 
pears due to the transition from the surface valence state 
So to the surface conduction state Ci. Considering the 
localization of the states (see Fig.|4]), namely the position 
of the adatom 3, where the charge density of So has been 
found, and that of the atoms 4-7-8 (the location of Ci) 
one sees that they are much nearer in y direction than 
in X. Therefore the transition probability in this case is 
higher for y polarized light, resulting in a negative RAS 
signal. Note that the transition from So to C2 is forbid- 
den. Also the peak P2 in the RAS is due to transitions 
between surface bands: the transitions from Si to both 
the surface conduction bands Ci and C2 are responsible 
for it. Here, the distance between adatom 9 (Si) and 
the filled pentagon 4-5-6-7-8 (Ci and C2) in a; direction 
is small, where also a little overlap of the charge den- 
sity appear. This yields in a strong positive signal in the 
RAS. Note that due to the overlap of Si with the bulk 
bands there is also a tiny contribution of bulk-surface 
transitions to P2. The peak P3 is determined by tran- 
sitions from S2 to Ci and C2. The vicinity of the state 
S2, which is located at atoms 1-2-3, to the states Ci and 
C2 at the pentagon with interstitial (4-5-6-7-8) is mainly 
due to its localization at adatom 3. From the charge 
density localized at the atoms 1-2, for both polarizations 
the transition probability should be similar since the dis- 
tances in X and y directions are equal. The additional 
localization at adatom 3 favor the y direction yielding 
the negative RAS signal. For the remaining two peaks 
P4 and P5 a mixture of bulk-surface and surface-bulk 
transitions (not restricted to single surface bands) is re- 
sponsible, therefore a clear picture could not be drawn. 
The peaks P2 and P3 are equal in height, whereas PI is 
a lower one. This might be due to the fact that for P2 
and P3 transitions to two surface conduction states are 
responsible, but one of them is forbidden for PI. How- 
ever, the overall magnitude of the probability of single 
transitions between surface bands is nearly the same in 
this case. 

Considering the dispersion of the eigenenergies, the re- 
gions in the IBZ with the main contributions to the RAS 
have been identified using again a single-k analysis^*': For 
PI, k points near F carry the main contribution, whereas 
for P2 k points close to the J-K boundary (large com- 
ponent), and for P3 those close to the F-J boundary 
(small ky component) comprise the main contributions. 
Such a trend has not been found neither for P4, nor for 
P5. A posteriori, due to this anisotropic distribution of 
the contribution, the convergence with respect to the k 
points had to be expected to be a delicate one. 

Changing the surface reconstruction certainly would 
affect the RAS. We can speculate that removing, e.g., 
the interstitial atom (atom 18) would change the local- 
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ization of the states Ci and C2 and which will have an 
effect to all three peaks discussed here. Adding an inter- 
stitial to the empty pentagon would affect mainly state 
S2, meaning the third peak in the spectra. Furthermore, 
having two nearly equivalent pentagons on the surface 
would influence the peak height, because the transition 
probability for the x and y would not be as significantly 
different as for the ADl reconstruction. However, also 
the energetic position of the surface peaks would indeed 
be affected. 



structure at around 3.2 eV and the double peak at 4.2 eV, 
which is shifted slightly to lower energies in our case. Also 
the experimental structures at around 5 eV can be found 
in our spectra, however, it is a little bit overestimated. 
The bulk derivative-like structure itself at around 3.5 eV 
is reproduced only qualitatively, which is a result of the 
small slab size. Compared to the tight-binding calcula- 
tion of Gavrilenko et al.^P, our spectra can resolve the 
characteristic low-energy surface peaks and describes the 
bulk derivative-like structure with the correct sign. 



C. Comparison with experimental results 

Finally, in Fig. [TO] we compare our RAS with the exper- 
imental data obtained by Aspnes and coworkers^^iM. Un- 
fortunately, the measured RAS did not show any surface- 
related spectral structure. This might be due to passiva- 
tion of the surface with Si-0 suboxide. Thus, the theo- 
retical surface peaks we have investigated cannot be com- 
pared. Regarding bulk-related structures, using a 11 DL 
slab does not allow the perfect description of the bulk 
derivative-like structure. In order to wash out artificial 
oscillations of bulk structures, we have applied a variable 
broadening. This means, we have employed a constant 
broadening for energies smaller than 2.5 eV (inset of the 
bulk spectra) and an additional broadening which goes 
linearly with the energy (factor 0.04) for higher energies. 
Furthermore, we scaled the measured RAS with a fac- 
tor of 10, which is the usually accepted factor for silicon 
surfaces In this way we got a reasonable agreement 
with respect to our slab thickness in the medium energy 
range. As visible, the calculation reproduces the spectral 
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FIG. 10: (Color online): Comparison of the RAS calculated 
using 25 k points in the 11 DL cell (light curves) obtained 
with a constant broadening of 0.08 eV (dotted line) and a 
variable one (solid line), see text. The experimental spectra 
are from Refs. [irlfi^ (dark dots connected by a solid line). A 
scissor operator of 0.5 eV has been applied to the theoretical 
spectra, and the RAS amplitude has been scaled by a factor 
of 10 (see text) . 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have investigated the electronic struc- 
ture and the RAS spectrum of the Si(113)3x2ADI. The 
electronic band structure and the DOS were found in 
excellent agreement with available experimental results. 
We identified three surface valence bands and two sur- 
face conduction bands near the Fermi level by computing 
the wave function localization in the surface region. The 
valence surface states are located at the adatoms and at 
the pentagon without an interstitial, whereas the conduc- 
tion surface states are located at the pentagon with an 
interstitial, in agreement with STM images. The three 
valence surface states in the energy gap compare with 
the two states which are visible in ARUPS"'^^ and the 
ARPES ^^'^^1^^ experiments, where in the latter one a 
third surface (resonance) state has been found in overlap 
with the bulk states. Nevertheless, the existence of three 
surface valence states has been confirmed by the compari- 
son of theoretical^i^^i^i^^ and experimental^i^i21i2aiSi2& 
STM images and by a core-level analysis.— Two of the 
valence surface states show a 3x2 and one a 3x1 peri- 
odicity, which has been confirmed experimentally. ^^'-^ In 
addition, we found a surface state below the valence band 
minimum. This state is located at the filled pentagon and 
shows a Tamm-like character. This finding calls for more 
experimental work for measuring the DOS at low energy. 

Concerning optical spectra, after checking the numeri- 
cal convergence with respect to the slab thickness and the 
number of k points, we have performed a layer-by-layer 
spectral decomposition. In this way the surface-relevant 
spectral features have been determined. Furthermore, an 
analysis of the contributions of the real and the imagi- 
nary part of the polarizability lead to the conclusion that 
the surface relevant structures are due to the imaginary 
part only. We have traced back the first three low-energy 
peaks of the RAS to transitions from the three valence 
surface states So, Si, and S2 (localized at the adatoms 
and the empty pentagon) to both the conduction sur- 
face states, Ci and C2 (localized at the pentagon with 
interstitial). These RAS peaks are hence very sensitive 
to a change of the surface reconstruction. The missing 
splitting of the peaks in the RAS with respect to the two 
surface conduction states is assumed to be due to the 
small energy difference between them. The only available 
experimental RAS does not show any low-energy surface- 



9 



related spectral features, probably due to a surface con- 
tamination. At higher energy, the bulk-derivative like 
structure, as well as other spectral features, agree with 
those of the measured RAS. Differences with respect to 
the experimental results can be due to the limited size of 
the slab used in the calculations and to the negligence of 
self-energy and excitonic effects. 
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